Self-trapping and stable localized modes in nonlinear photonic crystals 
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We predict the existence of stable nonlinear localized modes near the band edge of a two- 
dimensional reduced-symmetry photonic crystal with a Kerr nonlinearity. Employing the technique 
based on the Green function, we reveal a physical mechanism of the mode stabilization associated 
with the effective nonlinear dispersion and long-range interaction in the photonic crystals. 
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Photonic crystals are usually viewed as an optical ana- 
log of semiconductors that modify the properties of light 
similar to a microscopic atomic lattice that creates a 
semiconductor band-gap for electrons It is there- 
fore believed that by replacing relatively slow electrons 
with photons as the carriers of information, the speed and 
band-width of advanced communication systems will be 
dramatically increased, thus revolutionizing the telecom- 
munication industry. To employ the high-tech potential 
of photonic crystals, it is crucially important to achieve a 
dynamical tunability of their band gap . This idea can 
be realized by changing the light intensity in the so-called 
nonlinear photonic crystals, having a periodic modula- 
tion of the nonlinear refractive index |^ . Exploration of 
nonlinear properties of photonic band-gap (PBG) mate- 
rials may open new applications of photonic crystals for 
all-optical signal processing and switching, suggesting an 
effective way to create tunable band-gap structures op- 
erating entirely with light. 

One of the important ideas to control all-optical 
switching in the nonlinear regime is to explore the possi- 
bility of nonlinearity-induced self-trapping and nonlinear 
localized modes in photonic crystals. Existence of nonlin- 
ear localized modes for the frequencies in the forbidden 
gaps is usually associated with gap solitons, studied for 
one-dimensional Q and even two-dimensional (2D) mod- 
els Q , in the framework of the coupled-mode theory. Va- 
lidity of the coupled-mode theory is usually restricted by 
a weak modulation of the refractive index (the so-called 
shallow- grating case), and therefore this theory is not di- 
rectly applicable to the PBG crystals where modulation 
of the refractive index is of the order of its mean value. 
This observation calls for a systematic analysis of the 
(still open) problem of stable self-trapping and nonlinear 
localized modes in PBG materials, where the effects of 
discreteness j6| and long-range interaction have been 
recently shown to be of crucial importance. 

Nonlinear localized modes (also called intrinsic local- 
ized modes or discrete breathers) are associated with the 
energy localization that may occur in the absence of 
any disorder and solely due to nonlinearity |8j. Such 
nonlinear modes can be easily identified as approximate 



(or sometimes exact) analytical solutions of coupled- 
oscillator nonlinear lattice models and also in numeri- 
cal molecular-dynamics simulations, but only very re- 
cently the first observations of spatially localized non- 
linear modes have been reported in the physical systems 
of a very different nature ||^. The main purpose of this 
Letter is to predict the existence of nonlinear localized 
modes, analogous to gap solitons in the continuum limit, 
in 2D nonlinear photonic crystals, and to describe their 
unique properties including stability. 

We study nonlinear properties of 2D photonic crystals, 
assuming that their symmetry is reduced by inserting the 
rods made from a Kerr-type nonlinear dielectric material 
characterized by the third-order nonlinear susceptibility 
X^'^^. Specifically, we consider a periodic square lattice 
with the lattice spacing a which consists of two types 
of infinitely long cylindrical rods: the rods of radius ri 
made from a linear material and placed at the corners of 
the lattice, and the rods of radius r2 (r2 < ri) made from 
a nonlinear material and placed at the center of each unit 
cell (see top right inset in Fig. |l]). Linear properties of 
such photonic crystals are known [p^ . 

We assume that the rods are parallel to the X3 axis, 
so that the system is characterized by the dielectric con- 
stant e(x) — e(xi,X2)- In this case the evolution of the 
i?-polarized (with the electric field £J||a;3) light propa- 
gating in the {xi , a:2)-plane is governed by the scalar wave 
equation 



V'E{x,t)-\df[e{x)E]=0, 



(1) 



where V = d^_^ -\- d^^ and E is the x^ component 
of E. Taking the electric field in the form E{x,t) = 
g-iwt ]t;(x, t\u:) , where E{x, < | w) is a slowly varying en- 
velope, i.e. dfE{x, 1 1 oj) < ujdtE{x, t\uj), Eq. ([!]) can be 
reduced to 



E{x,t\ 



-2ie{x) 



Lo dE 
? 'dt 



(2) 



In the stationary case, i.e. when the r.h.s. vanishes, 
Eq. (^ reduces to the eigenvalue problem which can be 
solved in the linear limit i.e. when e{x) does not de- 
pend on the light intensity. In this case, the frequency 
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spectrum has a band-gap structure shown in Fig. It 
supports two band gaps, the lower of which extends from 
tu = 0.426 X 27rc/a to w = 0.453 x 2nc/a. 




FIG. 1. Band-gap structure of the _E-polarized light in 
a 2D reduced-symmetry photonic crystal with ri = 0.1a, 
r2 = 0.05a, and e = 11.4. Full lines are calculated directly 
from Eq. (Q) by MIT Photonic-Bands program whereas 
dashed line is found from linearized Eq. (5). The top right 
inset shows a cross-sectional view of the photonic crystal and 
the bottom right inset shows the corresponding Brillouin zone. 



A low-intensity light cannot propagate through a pho- 
tonic crystal if the light frequency falls into a band gap. 
However, it has been recently suggested [D that in the 
case of a 2D periodic medium with a Kerr-type nonlin- 
ear material, high-intensity light with the frequency in- 
side the photonic gap can propagate in the form of finite 
energy solitary waves — 2D gap solitons. Such solitary 
waves were analyzed in the framework of the coupled- 
mode continuum-limit equations valid for a weak modu- 
lation of the dielectric constant t{x). However, in real 
photonic crystals the modulation of e(a;) is comparable 
to its average value, so that the existence and stability of 
localized modes in such structures is still an open prob- 
lem. 




FIG. 2. Coupling coefficients Jn,m('^) for the photonic 
crystal depicted in Fig. |l|. The frequency uj — 0.4456 x 2-Kc/a 
is in the lower band gap. The contribution from Jo,o = 0.039 
is not shown. 



More specifically, the coupled-mode equations are valid 
if and only if the band gap Aw^ is vanishingly small, i. 
e. Acj^ ^ where A is an effective amplitude of the 
mode that is treated in the multi-scale asymptotic expan- 
sions as a small parameter. If we apply this model 
to describe nonlinear modes in a wider gap (see, e.g., 
discussions in Ref. we obtain a 2D cubic nonlinear 

Schrodinger (NLS) equation known to possess no stable 
localized solutions. Moreover, 2D localized modes of the 
coupled-mode equations are expected to possess an oscil- 
latory instability recently discovered for a broad class of 
ID coupled-mode equations ||l^]. Thus, it is clear that, if 
nonlinear localized modes do exist in realistic PEG ma- 
terials, their existence and stability should be associated 
with different physical mechanisms not accounted for by 
simplified continuum coupled- mode models. 

To study the nonlinear modes in such structures, we 
consider the nonlinear rods of small radius 7-2 as "defects" 
embedded into the linear photonic crystal formed by a 
square lattice of the rods of larger radius ri (" diatomic 
crystal"). Then, writing the dielectric constant e{x) as 
a sum of two periodic terms, €{x) = €i{x) -\- e2{x \ E) , 
where ei{x) describes the linear photonic crystal and 
£2(3; I E) corresponds to a lattice of nonlinear defect rods, 
one can present Eq. (||) in the form 

E{x,t\Lo)^ J d^y Gix,y\Lo)CE{y,t\Lo) , (3) 

where we introduce the Green function G{x, y | cj) of the 
linear photonic crystal (see, e.g., Ref. for its proper- 
ties) and the linear operator 

C^(^'^ye,{x\E) + 2^e{x)^^^. (4) 

Now, using the indices n and m for numbering the non- 
linear rods in the xi and X2 directions, we can describe 
their positions by the vectors Xn,m — nai+ma2, where 
ai and a2 are the primitive lattice vectors of the 2D pho- 
tonic crystal, and write 

e2{x I E) = {ef + x^''^\Eix, 1 1 ^ 0(a; - a;„,„) , 

n.Tn 

where 0{x) — 1 for \x\ < r2 and 0{x) — otherwise. The 
parameter e2*'' is the dielectric constant of the defect rods 
in the linear limit, while the term x'-^-'l^'P takes into ac- 
count a contribution due to the Kerr nonlinearity. If the 
radius r2 of the defect rods is sufficiently small, the elec- 
tric field E{x,t\LL!) inside them is almost constant, and 
Eq. (^ can be approximated by the discrete non- 
linear equation 

d 

Ot 

+ 5I'^"-fc.™-;M{4"^+x''^I^Ml'}^^fe,i = 0, (5) 

k,l 
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for the amplitudes of the electric field En^m{t\Lu) = 
E{xn,m,t\u) calculated at the defect rods. The param- 
eter a and the coupling coefficients 



(6) 



are determined by the Green function G{x,y\u!) which 
we calculate numerically by the FDTD method p!^ ] 
with the spatial step Ax = 0.01 to 0.03 and time step 
At = 0.005, for the lattice 24a x 24a. The Green func- 
tion and, therefore, the coupling coefficients Jn,m{^) 
in 2D photonic crystals are usually long-ranged func- 
tions. For instance, for the case of Fig. ^ we obtain 
J„,o - 0.012 (-1)" exp(-0.66 |n|) for n > 2, so that one 
should take into account the interaction between at least 
10 neighbors to achieve a good accuracy for the spectrum. 
By this means, Eq. (|^) is a nontrivial long-range gener- 
alization of the 2D discrete NLS equations extensively 
studied during the last decade ||l^ . We have checked the 
accuracy of the approximation provided by Eq. (H) by 
solving it in the linear limit. The low- frequency part of 
this dependence for 62*^^ = 11.4 is depicted in Fig. |l| by 
a dashed line; it has a minimum at w = 0.446 x 27rc/a 
being in a good agreement with the band edge calculated 
directly from Eq. (||). This lends a support to the valid- 
ity of Eq. (|^) and allows us to use this discrete model for 
studying nonlinear properties. 
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FIG. 3. (Color) Top and 3D views of a nonlinear local- 
ized mode in the 2D photonic crystal shown in Fig. |l| 
{uj = 0.442 X 2-Kc/a). 

Stationary nonlinear localized modes have been cal- 
culated numerically by the Newton-Raphson iteration 
scheme on the lattice 170a x 170a. Assuming that the lin- 
ear rods are made from GaAs whereas the nonlinear rods 



are made from some nonlinear material (which we do not 

specify varying e^2^), we find a continuous family of non- 
linear localized modes, and a typical example [smoothed 
by continuous optimization for Eq. (^] of such a mode 
is shown in Fig. ^. At first glance, this mode can be 
regarded as a donor state created by a single defect rod 
with larger dielectric constant. However, in the nonlin- 
ear case the mode stability becomes a critical issue. It 
can be determined from the so-called Vakhitov-Kolokolov 
stability criterion extended to the 2D discrete NLS mod- 
els in Ref. . According to this criterion, the nonlinear 
localized states with dQ/dio < are stable, and they are 
unstable otherwise. Here 



Q{uj) 



E 



(7) 



is the conserved mode power proportional to the energy 
of the electric field accumulated in the nonlinear localized 
mode. 
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FIG. 4. Power Q vs. frequency uj for the 2D nonlinear 
localized modes in the photonic crystal of Fig. 1 with two 
different . Solid lines - stable modes, dashed lines - un- 
stable modes. Grey areas show the projected band structure 
of the crystal. Insets show typical profiles of stable modes, 
and an enlarged part of the power dependence. 

As is well known (see, e.g., Refs. jl^) in the 2D discrete 
cubic NLS equation, only high-amplitude localized modes 
are stable, whereas no stable modes exist in the contin- 
uum limit. For our model, the high-amplitude modes are 
also stable (see inset in Fig. ^) , but they are not accessi- 
ble under realistic conditions: To excite such modes one 
should increase the refractive index at the mode center 
in more than 2 times. Thus, for realistic conditions and 
relatively small values of x^'^\ only low-amplitude local- 
ized modes become a subject of much interest since they 
can be excited in experiment. However, such modes in 
unbounded 2D NLS models are always unstable and ei- 
ther collapse or spread out ||l^. Here we reveal that, in a 
sharp contrast to the 2D discrete NLS models discussed 
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earlier in various applications, the low-amplitude local- 
ized modes of Eq. (|^) can be stabilized due to nonlinear 
long-range dispersion inherent to the photonic crystals. 
It should be emphasized that such stabilization does not 
occur in the models with only linear long-range disper- 
sion 

In order to gain a better insight into the stabiliza- 
tion mechanism, we have carried out the studies of Eq. 

for the exponentially decaying coupling coefficients 
Jn,m- Our results show that the most important factor 
which determines stability of the low-amplitude localized 
modes is a ratio of the coefficients at the local nonlin- 
earity (~ Ja,o) and the nonlinear dispersion (~ Jo.i). If 
the coupling coefficients Jn,m decrease with the distances 
n and m rapidly, the low-amplitude modes of Eq. (||) 

with e^°^ = 11.4 are essentially stable for Jo,o/'^o,i ^ 13. 
However, this estimation is usually lowered because the 
stabilization is favored by the presence of long-range in- 
teractions. 

It should be mentioned that the stabilization of low- 
amplitude 2D localized modes is not inherent to all types 
of nonlinear photonic crystals. On the contrary, the pho- 
tonic crystals must be carefully designed to support stable 
low- amplitude nonlinear modes. For example, in the pho- 
tonic crystal considered above such modes are stable at 
least for 11 < 62°'' < 12, however they become unsta- 
ble for e^"^ > 12 (see Fig. |). The stability of these 
modes can also be controlled by varying r2, ^i, or ei. 
Thus, experimental observation of the nonlinear localized 
modes would require not only the use of photonic materi- 
als with a relatively large nonlinear refractive index (such 
as GaAs/AlAs periodic structures or polymer PBG 
crystals [T9|]), but also a fine adjustment of the parame- 
ters of the photonic crystal. The latter can be achieved, 
in principle, by employing the surface coupling technique 
that is able to provide coupling to specific points of 
the dispersion curve, opening up a very straightforward 
way to access nonlinear effects. 

In conclusion, we have developed a consistent theory 
of nonlinearity- induced self-trapping effects in 2D nonlin- 
ear photonic crystals and predicted the possibility of the 
energy localization near the bandgap edge in the form of 
stable 2D nonlinear localized modes. 
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